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We propose a quantum sensing scheme for measuring weak forces based on a symmetry-breaking adiabatic 
transition in the quantum Rabi model. We show that the system described by the Rabi Hamiltonian can serve as 
a sensor for extremely weak forces with sensitivity beyond the yN /v/Hz range. We propose an implementation 
of this sensing protocol using a single trapped ion. A major advantage of our scheme is that the force detection 
is performed by projective measurement of the population of the spin states at the end of the transition, instead 
of the far slower phonon number measurement used hitherto. 
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Introduction. Using nanoscale mechanical oscillators as 
detectors of extremely weak forces has attracted considerable 
experimental interest [1]. Such systems allow to measure 
forces with sensitivity below the attonewton range which is 
beneficial for a broad range of applications. For example, a 
force detector with a nanomechanical oscillator coupled to a 
microwave cavity can reach sensitivity below one aN (10^^^ 
N) per v/Hz [2]. Other sensors use mechanical oscillators 
made of carbon nanotubes for force detection with sensitiv¬ 
ity in the zN (10^^' N) per y/Uz range [3]. Recently, the 
detection of ultra-weak forces as small as 5 yN (10^^^ N) 
was experimentally demonstrated using injection locked ions 
[4]. Force measurement in an ensemble of ions in a Penning 
trap uses Doppler velocimetry technique to detect force with 
sensitivity of 170 yN/y/Hz [5]. Another approach uses high- 
precision ion position determination to measure light pressure 
forces [6]. In all cases the force sensing based on mechan¬ 
ical oscillators is carried out by converting the force into a 
displacement that is measured by electrical or optical means. 

In this Letter, we introduce a different sensing protocol, 
which uses a system described by the quantum Rabi (QR) 
model as a probe that is sensitive to extremely weak forces. 
The QR model consists of a single bosonic mode and an ef¬ 
fective spin system which interact via dipolar coupling. We 
show that the effect of symmetry breaking of the underlying 
parity symmetry in the QR model due to the presence of exter¬ 
nal perturbation can be used in an efficient way for detection 
of classical oscillating forces. Our scheme relies on the adi¬ 
abatic evolution of the ground state of the QR model into the 
Schrodinger cat state, where the relevant force information is 
mapped in the respective probability amplitudes. The force 
sensing is performed simply by measuring the spin popula¬ 
tions. Therefore our protocol, which demands a single popu¬ 
lation measurement, is considerably faster than previous pro¬ 
tocols based on the detection of the motional degree of free¬ 
dom via Rabi oscillations. 

We consider a particular implementation of our sensing 
scheme using a coherently manipulated single trapped ion. 
The scheme, however, can be realized with various quantum 
optical systems such as nitrogen vacancy centers in diamond 
and superconducting qubits inside a microwave cavity [7-9]. 


We show that with the current ion trap technologies force sen¬ 
sitivity below one yN/v/Hz can be achieved. In addition, we 
show that our method can be applied for detection of spin- 
dependent forces which are created in magnetic-field gradi¬ 
ents or Stark-shift gradients. Hence our method can be used 
for studying magnetic dipole moments of atomic or molecular 
ions. 

Adiabatic quantum metrology using the quantum Rabi 
model. Our system consists of a two-level atom with states 
It) and It) coupled to a single bosonic mode described by the 
quantum Rabi model, 

HR = Ticod^d + ^^^^^^ay + Tigax{d^+ d). (1) 

Here and d are the creation and annihilation operators of 
bosonic excitation with frequency co and ap (j3 =x,y,z) are 
the respective Pauli matrices. The time-dependent Rabi fre¬ 
quency of the transverse field is given by £ly{t) and g is the 
spin-boson coupling. Recently, it was shown that the Rabi 
model permits exact integrability [10]. 

The quantum Rabi Hamiltonian (1) possesses a discrete 
symmetry revealed by the parity transformation d —d, 
O’, —O’, and (Jx-A- —Ox- In the following we consider the QR 
model in the ultrastrong coupling regime g/o < 1 and study 
the effect of a small perturbation term i/pert, which breaks the 
underlying parity symmetry of the model. The total Hamilto¬ 
nian including the perturbation term becomes 

H = HR + Hpert. (2) 

As we will see, the symmetry-breaking process allows us to 
estimate the perturbation term very accurately. 

The measurement protocol for i/pen starts by preparing the 
system in the ground-state of the unperturbed Hamiltonian Hr 
in the limit 12^(0) ^ g, such that |i//g(0)) = |—)y|0) where 

|-)y = (lt)-i|i))/V2and|«) is the Fock state of the bosonic 
mode with occupation number n. Then we adiabatically de¬ 
crease the transverse field fiyjf) in time such that the system 
evolves into the Schrodinger cat state 

Irg(O) =c+(f)|r+)+c_(f)|v/-), 


(3) 
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Figure 1: (Color online) (a) Low-energy spectrum of the quantum 
Rabi model as a function of time t. The two lowest lying states are 
separated by energy splitting Agap. At the initial moment ? = 0 the 
system is prepared in the ground-state with Fiy(O) 3> g and then it 
evolves into the Schrbdinger cat state (3). The energy difference be¬ 
tween the third excited energy (dashed line) and the ground-state en¬ 
ergy is Age. We assume y = 1.5 kHz, g = 25 kHz, to = 45 kHz and 
£2y(0) = 225 kHz. (b) Adiabatic parameter e versus time t for differ¬ 
ent values of y. 


where c±(f) are the respective probability amplitudes. Here 
11 /+) = |+)a: let) and I v/_) = |—)^ | —a) form the ground-state 
multiplet with |±)j. = (|t) ± \\))/\/2 and |a) stands for a co¬ 
herent state with amplitude a = —gjo). State (3) implies that 
for iTpert = 0 the parity symmetry is preserved by creating an 
entangled state with equal probabilities, c± = ±l/v/2. The 
effect of the perturbation is to break the parity symmetry of 
Hr by creating a ground-state wavefunction (3) with unequal 
probability amplitudes, |c+|^ ^ k-P- By measuring the re¬ 
spective probabilities at the end of the process one can esti¬ 
mate the unknown perturbation. 

In order to describe the creation of the symmetry-broken 
ground state we represent the Hamiltonian (2) within the 
ground-state multiplet, which leads to an effective two-level 
problem with the Hamiltonian [11] 


^eff = 


{W +1 -^pert I Y+) ^'^gap/2 

^ Agap /2 ( 1]/— I Hpeit I W— ) 


(4) 


Here Agap is the energy gap of the ground-state multiplet 
which scales with the transverse field as Agap ^ fiy (see. 
Fig. la) [11]. Hereafter we assume an exponential de¬ 
cay of the transverse field Cly{t) = ny( 0 )e^^ with a char¬ 
acteristic slope 7 which implies that Agap ~ The adia- 
baticity of the process is characterized by the condition e = 
I (v/g| Ii/e)/Age| ^ 1 , which requires the coupling be¬ 
tween the ground-state |t)/g) to the first excited-state jy/e) to 
be much smaller that the energy gap between them Age at any 
instant of time. In Fig. 1(b) we show the adiabatic parameter e 
during the creation of the Schrbdinger cat state (3) for various 
y. We observe that £ <C 1 and the non-adiabatic transitions are 
weaker for lower y. 

At the final stage the externally applied perturbation is de¬ 


tected by measuring the expectation value of < 7 ^. An analytical 
expression for the measured signal can be derived by solving 
the quantum evolution of a two-level system with the Hamil¬ 
tonian (4) for the specific time-dependence of D.y{t). 

Sensing weak forces and displacements. In the follow¬ 
ing we consider a harmonic oscillator represented by a single 
trapped ion with mass m confined in a Paul trap with an axial 
trap frequency ( 0 -. We assume that the effective spin system 
of the ion is implemented by two metastable atomic levels |t) 
and II) with a transition frequency Oo. We describe the small 
axial vibrations of the ion by the following motional Hamilto¬ 
nian 


H^ = hcod^d, z = zo{d^ + d), (5) 

where and d are the respective phonon creation and an¬ 
nihilation operators and zo = y/Ti/2m(0- is the spread of the 
oscillator ground-state wave function [ 12 ]. 

Electric field sensing. The ability to control the motional 
and internal states with high accuracy makes the trapped ions 
a formidable experimental tool for electric-field sensing [5, 
6 , 13]. In contrast to the conventional methods which rely 
only on the detection of the motional degrees of freedom [13, 
14], here the relevant information is transferred directly into 
the spin degrees of freedom due to the use of the symmetry¬ 
breaking adiabatic transition. In the following, we assume 
that a classical oscillating force F{t) = F(]Cos(a)(]f) with an 
amplitude — the parameter we wish to estimate — and 
frequency Wd = (0- — (0 shifted from the axial trap frequency 
CO, by a small detuning co (co- » CO) is applied to the ion. The 
action of the force is to displace the motional amplitude of the 
ion’s vibrational oscillator described by 

HF=F{t)z{t) = (6) 

where we have neglected the fast-rotating terms. In order to 
implement the spin-boson term in the quantum Rabi Hamilto¬ 
nian ( 2 ) we assume that the ion is simultaneously addressed 
by bichromatic laser fields in a Raman configuration with a 
wave vector difference Ak along the z direction, which in¬ 
duces a transition between the spin states via an auxiliary 
excited state. By setting the laser frequency beatnotes (Or = 
too — W, + 0 ) (Ob = tOo + ( 0 - — to close to the red- and blue- 
sideband transitions of the vibrational mode (O-, the result¬ 
ing Hamiltonian in the Lamb-Dicke limit (77 ^ 1) becomes 
[12, 15] 

Hs-b = (7) 


where g = Qrj is the spin-phonon coupling with £2 being the 
two-photon Rabi frequency and 77 stands for the Lamb-Dicke 
parameter. The transverse field in Eq. (2) can be created by 
driving the resonant carrier transition between the internal 
spin states using a microwave or radio-frequency field, which 
yields 


Hy(t) 


my(t) 


(k^it)ai+h.c) 


my(t) 


(Jy. 


( 8 ) 
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Figure 2: (Color online) (a) Time-evolution of the expectation value 
of Ox (solid lines) for minimal detectable force (11) and different 
y. The dashed lines correspond to the asymptotic solution (10a). We 
assume q single ^'*Mg+ trapped ion with an axial trap frequency cOr = 
6.3 MFIz. The other parameters are set to g = 25 klTz, co = 150 kHz, 
ny(0) = 225 kHz. (b) Signal-to-noise ratio versus We compare 
the numerical solution of the time-dependent Schrodinger equation 
with the Hamiltonian (2) (dots) with the asymptotic expression given 
by SNR = sinh(;rgzoF’(j/(liyft))) (solid lines). 



y(kHz) 

Figure 3: (Color online) The sensitivity tlforceV^ versus the slope 
y for various values of CO. We assume a single ^^Mg+ ion with an 
axial trap frequency cOr = 6.3 MHz. The interaction time is set to 
tf = 14y^*. The other parameters are g = 25 kHz, f2y(0) = 225 kHz. 
The exact solution (dots) is compared with the analytical expression 
^IforceV^ = (solid line) where F™" is given by Eq. (11). 


Here D.y{t) is the time-dependent Rabi frequency and we 
set the driving phase to 0 = —njl. In the interaction pic¬ 
ture rotating at the frequency (O the total Hamiltonian H = 
Hs-b + Hy + Hp is given by Eq. (2) where the symmetry¬ 
breaking term is 


Note that compared to other schemes here the sign of the 
force is fully preserved due to the tanh dependence of 
the signal. The corresponding signal-to-noise ratio SNR = 
((Jv(f/))/(A^(7j:(f/))^/^ is shown in Fig. 2(b). The minimum 
detectable force is determined by the condition that the signal- 
to-noise ratio is equal to one, which gives 


^pert = ^^(a’’'+fl)- (9) 

The force sensing protocol starts by initialization of the 
spins along the y direction and laser cooling of the single 
ion vibrational mode to the motional ground state. Subse¬ 
quently, the transverse held exponentially decays as Q.y{t) = 
nv(0)e^^, which drives the system adiabatically into the su¬ 
perposition state |v/g(f)) = c+(f) 11//+)-F c_(f) 11)/-). Here 
c±(f) are the respective probability amplitudes which are so¬ 
lutions of the time-dependent Schrodinger equation with the 
Hamiltonian (4) (see the Appendix). In our scheme the de¬ 
tection of the force is performed either by measuring the ex¬ 
pectation value of (Jx or by measuring the position quadrature 
Z = + d of the bosonic held. 

For vanishing force {F^ = 0) the parity symmetry is re¬ 
stored by creating an entangled ground state that is an equal 
superposition of states | V/±), which leads to {ax{tf)) = 0 and 
{Z{tf)) = 0. For yF 0, the parity symmetry of is bro¬ 
ken, which allow us to estimate by measuring OF or Z. In 
Fig. 2(a) the time-evolution of the expectation value of Gx 
in the presence of the symmetry-breaking term (9) is shown. 
At the interaction time tf ^ y the signal and its variance are 
described by (see Appendix) 

{Gx{tf)) = tanh (10a) 

{A^Gx{tf)) = l-{Gxitf))\ (10b) 


F7" = -^sinh-i(l). (11) 

ttgzo 

The sensitivity of the force measurement is dehned as T7foice = 
F™'^/y/v, where v = T/t is the repetition number with T be¬ 
ing the total experimental time. The time T includes the evo¬ 
lution time as well as the preparation and measurement times. 
Because our scheme relies on the projective measurement of 
the spin populations at the end of the adiabatic transition we 
have z Kitf. The sensitivity characterizes the minimal force 
difference, which can be discriminated within a total experi¬ 
mental time of one second. In Fig. 3 we show the sensitivity as 
a function of the slope y for various values of (0. Fowering y 
implies a longer interaction time tf and thus better sensitivity. 
For example, using the parameters in Fig. 3 for (£> — A5 kHz 
and interaction time tf = 30 ms we estimate force sensitivity 
of about 0.3 yN/v/Hz. Further increasing of the sensitivity 
to 0.16 yN/ v/Hz can be achieved with the interaction time of 
tf = 100 ms. 

Position quadrature. Alternatively, the force estimation can 
be carried out by measuring the expectation value of the posi¬ 
tion quadrature (Z(fy)). We hnd 




{A^Z{tf)) = \+A^-{Z{tf))\ 


( 0 ^ 


( 12 ) 
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Figure 4; (Color online) Signal-to-noise ratio versus the heating rate 
for various values of 7 and F^. We integrate numerically the master 
equation (14) with the Flamiltonian (2). The other parameters are set 
to g = 25 kHz, £0 = 150 kHz and f2v(0) = 225 kHz. 


Using Eq. (12) we obtain that for co > 2g the uncertainty of 
the position quadrature is higher than the measured signal 
(SNR < 1). At o = 2g and in the limit ^ 2yhl (ttzo), the 
SNR tends asymptotically to one from below and thus no mea¬ 
surement is possible, whereas for co <2g the force estimation 
is bounded by 


77 min 


Tiyco 

^gzo 


tanh 




(13) 


Comparing Eqs. (11) and (13) we conclude that the signal of 
(Jv provides better sensitivity. 

Effect of motional heating. The main source of decoherence 
which limits the force estimation in our scheme is caused by 
the motional heating. In order to account for the effect of 
the motional heating in the sensing protocol we numerically 
integrate the master equation 

^ = —-[H,p] + ^^^{h + l){2dpd^ — d^ap—pd^a) 
dt h 2 

+ ^^^h{2d^ pd — dd^ p — pdd"^). (14) 

Here /dec is the system decay rate and h is the average num¬ 
ber of quanta in the reservoir. In the limit « 1 the motional 

heating is characterized by the time f^ec = 1 /(K7dec) and the 
heating rate is (h) = 1 /fdec- In Fig- 4 we show the signal-to- 
noise ratio as a function of (h). By increasing the heating 
rate, the corresponding signal-to-noise ratio decreases with 
stronger damping for lower 7 and (O. Note that for small 7 
the effect of the motional heating is stronger, since the evolu¬ 
tion time f ~ 7^* is longer. Using the parameters presented 
in Eig. 4 we estimate sensitivity of 1.9 yN/v/Hz within the 
interaction time tf = 14 ms and heating rate approximately of 
0.1 ms^^ which corresponds to typical heating rates in linear 
ion Paul traps. Eor a cryogenic ion trap with heating rate of 
order of (n) = 0.01 ms^* [16] the force sensitivity would be 


approximately of 0.7 yN/v/Hz for co = 15 kHz and 7 = 0.8 
kHz. 

Sensing of spin-dependent forces. Here we introduce a 
measurement protocol that can be used for sensing of spin- 
dependent forces. The origin of such forces can be due to 
(i) spatially dependent Zeeman shifts experiences by spins 
in a magnetic-field gradient [17] or (ii) spatially dependent 
Stark shifts experienced by atoms in an intensity field gradi¬ 
ent, which give rise to an optical dipole force [18]. Consider 
the case of an oscillating magnetic field gradient along the 
trap axis. Consider also an ion chain with two ions where 
the internal states of the first ion are formed by the magnetic 
insensitive clock states with s =t,i- These clock states 
are used to implement the quantum Rabi Hamiltonian (1) tak¬ 
ing Ox = (It)ci(il +h-c-) and C7v = i(|i)ci (t| - h.c.) with the 
same method as was described above. The auxiliary ion (not 
necessarily the same atomic species) is prepared in one of its 
magnetic sensitive states |aux). The external oscillating mag¬ 
netic field gradient B{t) = cos{cot)Bozez will create a coupling 
between the auxiliary ion states and the collective vibrational 
modes. Assuming that the oscillating frequency is close to the 
center-of-mass vibrational mode (0 = (Oc.m. — (0 and neglect¬ 
ing the fast oscillating terms we arrive to 


Tfpert — 


Zc.m.Fi 

~7T' 


(fll -F fl)|aux)(aux|, 


(15) 


where Zc.m. = \/fi/2mcUc,m.- Here 77 = gjp^BQ/2 is the force 
acting on the magnetic moment = (.g7/iB/2)|aux)(aux| 
of the auxiliary ions where .gy is the Lande g-factor and pg is 
the Bohr magneton. The force acting on the clock ion depends 
on the spin state of the auxiliary ion such that the perturbation 
term becomes 

(aux|Hpert |aux) = (ig) 

v2 

The induced force on the clock ion can be used for studying, 
e.g., the magnetic moment of the auxiliary atomic or molecu¬ 
lar ions trapped together with the clock ion. Indeed, assuming 
magnetic dipole moment 2pg of the auxiliary ion and mag¬ 
netic field gradient of 1 T/m the resulting force is about 9 yN, 
which can be detected by measuring the spin population of the 
clock ions. 

Conclusions. We have shown that the system described 
by the quantum Rabi Hamiltonian can serve as a detector of 
extremely small forces. The underlying physical mechanism 
is the process of symmetry-breaking adiabatic transition due 
to the presence of force perturbations. Our sensing protocols 
can be implemented using a trapped ion, where the parameters 
which drive the system across the adiabatic transition are con¬ 
trolled by external laser or microwave fields. We have shown 
that a system of a single trapped ion can be used as a probe 
for electric sensing with sensitivity about and even below 1 
yN/v/Hz range. Additionally, the proposed method can be 
extended for sensing magnetic fields. A major advantage of 
our protocol is that it demands a single population measure- 
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ment, thereby achieving a considerably speed-up over previ¬ 
ous protocols using phonon number measurement via Rabi os¬ 
cillations. 

This work has been supported by the EC Seventh Frame¬ 
work Programme under Grant Agreement No. 270843 (iQlT). 


APPENDIX: DERIVATION OE THE SIGNAL AND THE 
VARIANCE OE THE SIGNAL 


In the symmetry-broken phase the underlying dynamics of 
the QM Hamiltonian with the perturbation term (9) can be 
captured within the two-level model with the effective Hamil¬ 
tonian 


zoFii ^ /iAgap {t)/2 
^Agap(f)/2 


(17) 


The probability amplitudes c±{t) for the state vector 
\Vg{t)) = c+(f)IV^+) +^-(01*/-) obey the time-dependent 
Schrodinger equation, which turns into a system of two cou¬ 
pled differential equations, 

ic+(0 = -^c+(0 + ^^c_(f), (18a) 

= + (18b) 

subject to the initial conditions c+(f,) = l/-\/2 and = 

-1/V2. We assume that the transverse field varies in time 
as ily{t) = , which implies that Agap(f) = 

and the two-state problem reduces to the Demkov model. The 
probabilities at ff ^ 7^^ are given by 


\c+{tf)? 


1 

2 


■ tanh 


\ tiyco ) 


(19) 


and |c_(f/)p = 1 — |c+(f/-)p. The expectation value and the 
variance of Ox with respect to | ^re 

{ax{tf))=2\c+{tf)\^-l, (20a) 

{A^aitf))=4\c+{tf)\^\c-{tf)\\ (20b) 

For the measured signal of the quadrature Z and its variance 
we find 
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